Analogues of the adjoint matrix for generalized 
inverses and corresponding Cramer rules. 
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Abstract 

In this article, we introduce determinantal representations of the 
Moore - Penrose inverse and the Drazin inverse which are based on 
analogues of the classical adjoint matrix. Using the obtained analogues 
of the adjoint matrix, we get Cramer rules for the least squares solution 
and for the Drazin inverse solution of singular linear systems. Finally, 
determinantal expressions for A^A, AA+, and A^A are presented. 
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1 Introduction 

Determinantal representation of the Moore - Penrose inverse was studied in 
[H [21 [HI [m [I2]- The main result consists in the following theorem. 

Theorem 1.1 The Moore - Penrose inverse A"*" = (a^)„xm o/ A G CJ?^" 
has the following determinantal representation 



a, 



E 

(a,l3)eNr{j,i} 



a ^ 

daj i\ /3\ 



A' 



E 



I A"^! 
1^5 1 



1 < i, j < n. 
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Stanimirovic' [13] introduced a determinant al representation of the Drazin 
inverse by the following theorem. 

Theorem 1.2 The Drazin inverse = (a^) of an arbitrary matrix A G 
k possesses the following determinantal representation 



C"^" with IndA 



n 



E 

{a,^)GNr^{j,i} 



a 



E 

where s > k and = rankA'^. 



:a^ 



1^5 I 



l<i,j <n; 



These determinantal representations of generalized inverses are based on cor- 
responding full-rank representations. 

We use the following notations from [H [12]. Let C™^" be the set of m by 
n matrices with complex entries, C^^" be the subset of C™"^" in which every 
matrix has rank r. Im denotes the identity matrix of order m, and ||.|| = ||.||2 
is the Euchdean vector norm. Let a := {Q!i,...,afc} C {l,...,m} and 
/3 := . . . , (3k} C {1, . . . , n} be subsets of the order 1 < A; < min {m, n}. 
Then lAol denotes the minor of A G C™^" determined by the rows in- 



dexed by a and the columns indexed by (3. Clearly, 



denotes a prin- 



cipal minor determined by the rows and columns indexed by a. The co- 



factor of a 



is denoted by 



d 



For 1 < < T2,, de- 



in A e C 

note by Lfc,„ := {a : a = (ai, . . . , au) , 1 < cti < . . . < ak < n} the collec- 
tion of strictly increasing sequences of k integers chosen from {1, . . . , n}. Let 

Nk ■■= Lk,m X Lk-n- For fixed a G Lp^m, P G Lp^n, 1 < p < k, let 

h,m (a) ■= U ■ I & Lk,m, I ^ a}, 
■.= {J : J e Lfc,„, J D /?}, 
Nk := Ik,m {a) x Jk,n iP) 

For case i E a and j E (3, we denote 

h,m{i} := {a : a G Lk,m,i e a},Jk,n{j} ■= {(3 : [3 e Lk,n,j e /?}, 

Nk{i,j} ■= h,m{i} X Jk,n{j}- 

In this paper we introduce determinantal representations of the Moore - 
Penrose inverse and of the Drazin inverse based on corresponding limit rep- 
resentations. The obtained determinantal representations can be considered 
as founded on some analogues of the classical adjoint matrix. The corre- 
sponding Cramer rules for the complex system of linear equations with a 
rectangular or singular coefficient matrix follow from these analogues. 



2 



2 Analogues of the classical adjoint matrix 
for the Moore - Penrose inverse 

We shall use the following well-known facts (see, for example, [9]). 

Definition 2.1 The matrix G C"^™" is called the Moore - Penrose in- 
verse of an arbitrary A G C™^" if it satisfies the equations 

AA+A = A; A+AA+ = A+; (AA+)* = AA+; (A+A)* = A+A. 

The superscript * denotes conjugate transpose matrix. 

Lemma 2.1 ^ There exists a unique Moore - Penrose inverse A+ o/ A G 

Lemma 2.2 ;9] If A e C™^", then 

A+ = limA* (AA* + AI)~^ = lim (A*A + AI)"^ A*, 

A^O A-i-O 

where A G M+, and M_|_ is the set of positive real numbers. 

Lemma 2.3 /5/ If A E C™^", then the following statements are true, 
i) If Tank A = n, then A+ = (A*A)"^ A* . 
a) //rank A = m, then A+ = A* (AA*)"^ . 
Hi) If Tank A = n = m, then A+ = A^"*^ . 

Theorem 2.1 Let dr be the sum of principal minors of order r o/ A G 
C"^". Then its characteristic polynomial p a {t) can be expressed as pa (t) = 
det {tl - A) = r - dit"-! + £^2^-2 - . . . + (-1)" dn. 

Denote by a j and cij the jth column and the ith row of A G C"^" respec- 
tively. In the same way, denote by a*- and a* the jth column and the ith row 
of Hermitian adjoint matrix A*. Let A.^ (b) denote the matrix obtained from 
A by replacing its jth column with some vector b, and let Aj, (b) denote the 
matrix obtained from A by replacing its ith row with b. 

Lemma 2.4 If A e C™^", then Tank (A*A) . (a*^.) < r. 
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Proof. Let Pjfc(— a^fe) G C"^", {k ^ i), be the matrix with —ajk in the 
{i, k) entry, 1 in all diagonal entries, and in others. It is the matrix of an 
elementary transformation. It follows that 



( E otfeOfei • • • «y ■ ■ ■ E a*lk<^kn \ 



E «nfe«fel • • • «nj ■ ■ ■ E Onfe^fen 
i—th 



The obtained above matrix has the following factorization. 



/ E «ife«fei 
E <it«fei 



E «lfc«fcn \ 

ky^j 



i—th 



a 



12 



Ooi a. 



21 "22 • • 
\ "^nl ^n2 • • 

/ an 



Denote by A := 







"-Im 
«2m 

^nm / 




an 





V Oml 



E anfe^fc" 
k^j 





1 





i—th 











-th 



j — th. The matrix A is ob- 



tained from A by replacing all entries of the jth row and of the ith column 
with zeroes except that the (j, i) entry equals 1. Elementary transformations 
of a matrix do not change its rank. It follows that rank (A* A) ^ (a* ) < 

min |rank A*, rank a|. Since rank A > rank A = rank A* and rank A* A = 
rank A the proof is completed. ■ 

The following lemma is proved in the same way. 



Lemma 2.5 If A e C^"^", then rank (AA*) . (a*.) < r. 
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Theorem 2.2 The Moore-Penrose inverse A+ of A E CJ!*^" can be repre- 
sented as follows 



h 



^dr(A*A) 

E 

/3eJr, n{«} 



where 



k,= E,, ((A*A).,(a*.))^ rf.(A*A)= E (A* A) 



or 



A^ 



, where 



dr{AA*) 

E |((AA*),.(a*,)):|, 4(AA*)= E l(AA*) 



(2) 



a€lr,m{j} 



Proof. At first we shall obtain the representation ([T]). If A G then the 
matrix (AI + A* A) G C"^" is Hermitian and rank (AI + A* A) = n. Hence, 
there exists its inverse 



(AI + A*A) 



-1 



1 



/ Lii L21 
L12 L22 



det (AI + A* A) 



Lnl \ 
Ln2 



\ Lin L2n • • • Ln « / 



where Lij = l,n) is a cofactor in AI + A* A. By Lemma 12.21 A^ 

lim (AI + A*A)~^ A*, so that 



A+ = lim 

A-)-0 



/ det(AI+A*A) ;^(a*J 
dct(AI+A*A) 



det(AI+A'A)^(a'j; 
\ det(AI+A*A) 



dctCAI+A'A) i(a*^) \ 
det(AI+A*A) 



det(AI+A*A) „(a*^) 
det(AI+A*A) / 



(3) 



From Theorem 12.11 we get 

det (AI + A* A) = A" + rfiA"-^ + rfaA"-^ + . . . + d„ 



where dj. (Vr = 1, n — 1) is a sum of principal minors of A* A of order r and 
dn = det A* A. Since rank A* A = rank A = r, then dn = dn-i = ... = 
dr+i = and 



det (AI + A* A) = A" + diA"-' + dsA""" + . . . + dr^ 



(4) 
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In the same way, we have for arbitrary 1 < i < n and 1 < j < m from 
Theorem 12.11 



det (AI + A*A)^ (a*.) = /f^U"-^ + /J^^'^A^-^ + . . . + l^^^ 



where for an arbitrary 1 < < ri — 1, / 



= det (A*A) . (a*j.). By Lemma [231 rank (A* A)^ (a*^.) < r so that if 



E ((A*A).,(a*.)) 



3'l f} 



, and 



k > r, then 



((A*A).,(a*.))^ = 0, (WP e Jfe,nW,Vz = l,n,Vj = l,m). 

= and 



Therefore if r + 1 < < n, then 



4*^) = det (A*A) . (a*^.) = 0, (Vi = l,n, Vj = l,m). Finally we obtain 



E ((A* A)., (a*,.)) 



det (AI + A*A)^ (a*.) = /f ^'^A""^ + 4''^A"-2 + . . . + /fe)A 



(5) 



By replacing the denominators and the numerators of the fractions in 
entries of matrix (E]) with the expressions and (jS]) respectively, we get 



zJ">A"-i+...+«5.^^'A"-'' 
A"+diA"-i+...+drA"-'" 



(11). 



A^ 



lim 

A^O 



A" 



A"+diA"-i+...+drA"-'' 



A"-i+...+;("i)A 



* A"+(iiA" 



-...+drA" 

dr 



,(nl) 
, £r 

\ dr 



A"+diA"-i+...+(irA"-'' 



,(lm) 



(nm) 



From here the representation ([T]) of A+ follows by denoting Ir'''^ = I 
We obtain the representation ([2]) in the same way. ■ 



Remark 2.1 //rank A = n, then from Lemma WT^ we get = (A* A) ^ A*. 
Representing (A*A)"-'^ by the classical adjoint matrix, we have 



A^ 



^ /det(A*A).i (a*i) ... det(A*A).i (a*, 



det(A*A) \^det(A*A).„(a*i) ... det(A*A).„ (a*. 
If n <m, then ^ is valid. 



(6) 
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Remark 2.2 As above, z/ rank A = m, then 



A' 



1 



'det(AA*)i. (a* 



det(AA*) \^det(AA*)i. (a; 

If n> m, then ^ is valid as well. 



det(AA* 
det(AA* 



(7) 



Remark 2.3 The representation (QP can he obtained from Theorem \1.1\ by 
using the Binet-Cauchy formula. We use another method, which is the same 
for the determinantal representations of the Moore-Penrose inverse by ^ 
and and of the Drazin inverse by ( fTTl) . 



Remark 2.4 To obtain an entry of by Theorem M . 1\ one calculates (C^C^+ 
^n-i^m-i) determinants of order r. Whereas by (CP we calculate as much as 
(C^ + Cl^) determinants of order r or we calculate the total of (C^ + C^l\) 
determinants by Ij^. Therefore the calculation of entries o/A+ by Theorem 
\2.2\ is easier than by Theorem 



Corollary 2.1 If A E C™^" and r < min{m,n} or r 
projection matrix P = A+A can be represented as 

Pij 



m < n, then the 



dr{A*A)^ 

where d.j denotes the jth column of (A* A) and, for arbitrary 1 < i, j < n, 
p,,= E ((A*A).,(d,))^ . 

l3£Jr,n{i} 

Proof. Representing the Moore - Penrose inverse A"*" by 

/ill h2 ■■■ hm\ ( 0-11 ai2 

/21 I22 ■ ■ ■ hm 



1 

dr (A* A) 



\^nl ^n2 • • • ^nm/ 

Therefore, for arbitrary 1 < i, j < n we get 



^21 CL22 
\0-ml 0,m2 



we obtain 
ain \ 

02n 



P 



E E 

k fSeJr.n{i} 



((A*A).,(a*,))^ 



E E ((A*A).,(a:,-a,,))^ 



^eJr,n{*} k 



E ((A*A).,(d*.))^ 



/3eJr,n{*} 



Using the representation ([2]) of the Moore - Penrose inverse the following 
corollary can be proved in the same way. 
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Corollary 2.2 If A E CJJ^^", where r < min {m,n} or r = n < m, then a 
projection matrix Q = AA+ can be represented as 



Q 



dr (AA*] 



where gi denotes the ith row of (AA*) and, for arbitrary I < i, j < m, 
5.,= E |((AA*),,(gO)^|. 

a£lr,m{j} 

Remark 2.5 By definition of the classical adjoint Adj{A) for an arbitrary 
invertible matrix A E C*^^" one may put Adj (A) ■ A = det A-I„. If A G C"^" 
and rank A = n, then by Lemma \2. 3\ A+A = I„. Representing the matrix 
A"*" by ^ as A+ = det(A*A) ' '^^ obtain LA = det (A*A)-I„. This means that 
the matrix is a left analogue of Adj {A) , where A G C™^". //rank A = m, 
then by Lemma \2.3[ AA+ = I^- Representing the matrix A"*" by ^ as 
A+ = ^^g^j-^A*) ' obtain AR = • det (AA*). This means that the matrix 
R zs a right analogue of Adj{A), where A G C™^" . 

Remark 2.6 If A E C™^" and r < min{m,n}, then by (C]) we have A+ = 
^^^^j^, whereL = {kj) G C"^™. From Corollary\K^we getLA = 4 (A*A)- 
P. The matrix P is idempotent. All eigenvalues of an idempotent matrix 
chose from 1 or only. Thus, there exists an unitary matrix U such that 
LA = dr (A* A) Udiag (1, . . . , 1, 0, . . . , 0)U*, where diag (1, . . . , 1, 0, . . . , 0) G 
£_nxn ^ diagonal matrix. Therefore, the matrix L can be considered as a 
left analogue of Adj {A), where A G C^^". 

In the same way, if A E CJ!*^" and r < min{m,n}, then by ^ we 
have A+ = -j-jxa*)' '^^^^^ = {fij) ^ C"^™. From Corollary \2.2\ we get 
AR = (AA*) • Q. The matrix Q is idempotent. There exists an uni- 
tary matrix V such that AR = dr (AA*) Vdiag (1, . . . , 1, 0, . . . , 0)V*, where 
diag (1, . . . , 1, 0, . . . , 0) G C"^^™. Therefore, the matrix R can be considered 
as a right analogue of Adj (A) in this case. 

3 An analogue of the classical adjoint matrix 
for the Drazin inverse 

Definition 3.1 ^ \^ Let A G C"^" with IndA = k, where a nonnega- 
tive integer IndA := min |rankA'^+-'^ = rankA'^). Then the matrix X 

keNU{0} ^ ^ 
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satisfying 

A*-'+^X = A^ XAX = X; AX = XA (8) 

is called the Drazin inverse of A and is denoted by ^ = . In particular, 
if IndA = 1, then the matrix X in ^ is called the group inverse and is 
denoted by X. = A*. 

Remark 3.1 If IndA = 0, then A is nonsingular, and A^ = A^^. 

The Drazin inverse can be represented explicitly by the Jordan canonical 
form as follows. 

Theorem 3.1 gj/ // A G C"^" with IndA = k and 

where C is nonsingular, rankC = rankA^, and N is nilpotent of order k, 
then 




We use the following theorem about the limit representation of the Drazin 
inverse. 

Theorem 3.2 ^ // A G C"^", then 

A^ = lim (AI„ + A'^+i)"^ A^ 
where k = IndA and X G M+. 

Denote by a'-^'^'' and a^-'^'' the jth column and the ith row of A'^ respectively. 
Lemma 3.1 If A e C"^" with IndA = k, then 

rank A^+^ (^af^ < rank A^+\ ^i,j = hn. (10) 

Proof. The proof of this lemma is similar to that of Lemma [2.41 

In the following theorem we introduce a determinantal representation of 
the Drazin inverse. 
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Theorem 3.3 // IndA = k and rank A'^"'"^ = rank A'^ = r < n for an 

arbitrary matrix A G C"^", then 



D 



where 



(11) 



fSeJr.n{i} 



(Vz,j = l,n) 



Proof. The proof of this theorem is analogous to that of Theorem 12.21 by 
using Theorem I2.H Lemma I3.H and Theorem 13.21 

In the foUowing coroUaries we introduce determinantal representations of 
the group inverse A* and the matrix A^A respectively. 

Corollary 3.1 If IndA = 1 and rank A^ = rank A = r < n for an arbitrary 
matrix A G C"^", then 



A* 



9ij 



dr (A^] 



where 4 (A^) = ^ i^X 



l3€Jr,n{i} 



Proof. The proof follows from Theorem 13.31 in view of A; = 1. 

Corollary 3.2 // IndA = k and rankA'^^-'^ = rankA'^ = r < n for an 
arbitrary matrix A G C"^", then 



A^A 



dr f A'^+n 



where Vij 



E {Ar{^/'^'%,{^^,J = l;^). 



l3e.Jr,n{i} 

Proof. Representing the Drazin inverse A^ by flTTl) we obtain 

di 



A^A 



dr fA^-+ll 



dr{^'^')Jr.xn 
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Here for arbitrary 1 < i, j < n we have 



Vij = E E 

S /3eJr,n{«} 



E E 

/3eJr,n{j} S 



A 



a.s' ■ a 



E 



Remark 3.2 The matrix (A^A) is idempotent. Similarly to the case of 
Remark \2.6[ the matrix D = (dij)^^^ can be considered as an analogue of the 
classical adjoint matrix, where D = dr (A'^+-'^)A^. 



4 Cramer rules for generalized inverse solu- 
tions 



Definition 4.1 Suppose in a complex system of linear equations: 

A-x = y 



(12) 



the coefficient matrix A G C™^" and a column of constants y = {yi, . . . , ymf' G 
C™. The least squares solution of the system (W\) is the vector x° G C" sat- 
isfying 

||x°|l = min< ||x|| I IIA ■ x — yll = minllA ■ x — yll >, 

where C" is an n-dimension complex vector space. 

Theorem 4.1 The vector x = A+y is the least squares solution of the 
system H^) . 



Theorem 4.2 The following statements are true for the system of linear 
equations / f7^) . 

i) //rank A = n, then the components of the least squares solution x = 
(x'j', . . . , are obtained by the formula 



X 



_ det(A-A)., (f) 
det A* A 



(Vj = l,n), 



(13) 



where f = A*y. 
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a) //rank A = r < m < n, then 



E ((A*A).,(f))^ 
x° = '''-^^^ , (Vj = M). (14) 



dr (A* A) 



Proof, i) If rank A = n, then we can represent A+ by ([7]). By multiplying 
A+ into y we get f|T3|) . 

ii) If rank A = k < m < n, then A"*" can be represented by ([T]). By 
multiplying A+ into y the least squares solution of the linear system f[T^ is 
given by components as in f[T^ . ■ 

Using ([3]) and ([8]), we can obtain another representation of the Cramer 
rule for the least squares solution of a linear system. 

Theorem 4.3 The following statements are true for a system of linear equa- 
tions written in the form x ■ A = y. 

i) //rank A = m, then the components of the least squares solution x° = 
yA+ are obtained by the formula 

det(AA*),.(g) 

m 5 Vz = l,m , 

det AA* ^ ^ 

where g = yA*. 

ii) //rank A = r < n < m, then 

E l((AA*),.(g)):i 

aG-fr,m{«} / X 

Xi = , , . ■ , , I V« = 1, ml . 

' dr{AA*) ^ ' 

Proof. The proof of this theorem is analogous to that of Theorem 14.21 

Remark 4.1 The obtained formulas of the Cramer rule for the least squares 
solution differ from similar formulas m |3, 0, [721 [121 HE/- They give a 
closer approximation to the Cramer rule for consistent nonsingular systems 
of linear equations. 

In some situations, however, people pay more attention to the Drazin 
inverse solution of singular linear systems [6l [16] . Consider a general system 
of linear equations f[T2l) . where A e C"'^'^ and x, y are vectors in C". R{A) 
denotes the range of A and A^(A) denotes the null space of A. The charac- 
teristic of the Drazin inverse solution A^y is given in [16] by the following 
theorem. 
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Theorem 4.4 Let A G C^^^" with Ind{A) = k. Then A^y is both the 
unique solution in i?(A^) of 



A^-'+^x = AV, 

and the unique minimal P-norm least squares solution of ( fl^ 



(15) 



Remark 4.2 The P-norm is defined as ||x||p = ||P ^x|| for x G C", where 
P is a nonsingular matrix that transforms A into its Jordan canonical form 



Remark 4.3 Since l\T5\) is analogous to the normal system A* Ax = A*y, 
the system l[T5\) is called the generalized normal equations of [W^) . (see JT^). 

We obtain the Cramer rule for the P-norm least squares solution of ( 1T2|) in 
the following theorem. 

Theorem 4.5 Let A G C"^" with Ind{A) = k. Then the unique minimal 
P-norm least squares solution x = {xi, . . . of the system / fl^) is given 

by 

l3eJr, n{«} 



Xi 



Vi = 1, n, 



(16) 



where g = A'^y. 



Proof. Representing the Drazin inverse by ([TTj) and by virtue of Theorem 
14.41 we have 



Xi 



X 



A^y 



dr (A'^+l) 



[ E disys\ 

s=l 

n 

^E d-nsVs J 



Therefore, 



Xi 



1 " I / 
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/3eJr,n{i} 

From this f|T6l) follows immediately. 



5 Examples 

1. Let us consider the system of linear equations. 

2xi — 5x3 + = 1, 
7xi — 4x2 — 9x3 + 1.5x4 = 
3xi — 4x2 + 7x3 ~ 6.5x4 = 
Xi — 4x2 + 12x3 ~ 10.5x4 ■ 

The coefficient matrix of the system is the matrix A = 

Vl -4 

We calculate the rank of A which is equal to 3, and we have 



1. 



12 



(17) 



/2 





-5 


4 \ 


7 


-4 


-9 


1.5 


3 


-4 


7 


-6.5 



-10.5/ 



A* 



/ 2 


7 


3 


1 \ 





-4 


-4 


-4 


-5 


-9 


7 


12 


V4 


1.5 


-6.5 


-10.5/ 



, A*A 





63 


-44 


-40 


-11.5\ 




-44 


48 


-40 


62 




-40 


-40 


299 


-205 


V 


-11.5 


62 


-205 


170.75/ 



At first we obtain entries of A+ by (Q: 



dsiA*A) 



63 
-40 
-11.5 



63 -44 
-44 48 
-40 -40 
-40 -11.5 
299 -205 
-205 170.75 



-40 
-40 
299 



48 
-40 
62 



63 
-44 
-11.5 
-40 



-44 
48 
62 
62 



-11.5 

62 
170.75 



299 -205 
-205 170.75 



102060, 





2 


-44 


-40 




2 


-44 


-11.5 




2 


-40 


-11.5 


/ll = 





48 


-40 


+ 





48 


62 


+ 


-5 


299 


-205 




-5 


-40 


299 




4 


62 


170.75 




4 


-205 


170.75 



25779, 



14 



and so forth. Continuing in the same way, we get 



A^ 



102060 



/ 25779 
-3840 
28350 

y 39558 



-4905 20742 

-2880 -4800 

-17010 22680 

-18810 26484 



-5037 \ 
-960 
-5670 
-13074/ 



Now we obtain the least squares solution of the system ([T7|) by the matrix 
method. 



/ X 

X 

x' 



11 \ 


21 



V^41, 



102060 



/ 25779 
-3840 
28350 

\ 39558 



-4905 20742 

-2880 -4800 

-17010 22680 

-18810 26484 



-5037 \ 
-960 
-5670 
-13074/ 



/1\ 

2 

3 
VV 



102060 



/ 73158 \ 
-24960 
56700 
\ 68316 / 



/ 12193 \ 

1071 
5 

9 

\ 5693 ; 
\ 8505 / 



Next we get the least squares solution of the system (H 
([HI), where 



by the Cramer rule 



/2 


7 


3 


1 \ 






/26\ 





-4 


-4 


-4 


2 




-24 


-5 


-9 


7 


12 


3 




10 


U 


1.5 


-6.5 


-10.5/ 


v) 







Thus we have 



102060 



+ 



102060 



63 
-44 
-40 



26 


-44 -40 




26 


-24 


48 -40 


+ 


-24 


10 


-40 299 




-23 


-40 


-11.5 




73158 


299 


-205 


) 






-205 


170.75 




102060 



26 -40 
-24 -40 
10 299 



-44 -11.5 
48 62 
62 170.75 

12193 



17010' 



63 
-44 
-11.5 



26 
-24 
-23 



-11.5 

62 
170.75 



+ 



15 



-24 -40 62 
10 299 -205 
-23 -205 170.75 



102060 



63 
-44 
-40 



-44 26 
48 -24 
-40 10 



+ 



+ 



48 
-40 
62 



-24 
10 
-23 



102060 



63 
-44 
-11.5 



62 
-205 
170.75 

-44 26 
48 -24 
62 -23 



+ 



48 -40 -24 
-40 299 10 
62 -205 -23 





24960 


416 


102060 


1071' 




63 


26 -11.5 




-40 


10 -205 




-11.5 - 


-23 170.75 




56700 


5 




102060 


~ 9' 




63 


-40 26 


+ 


-40 


299 10 




-11.5 


-205 -23 


68316 


5693 


102060 


8505' 



2. Let us consider the following system of linear equations. 

Xi - X2 + X3 + X4 = 1, 

X2 - X3 + X4 = 2, 
xi — X2 + X3 + 2x4 = 3, 

Xl - X2 + X3 + X4 = 1. 



The coefficient matrix of the system is the matrix A 



It is easy to verify the following: 
/3 -4 4 3\ 



/I 


1 





4 

V3 



-5 
-4 



5 
4 





4 

3/ 



and rank A = 3, rank A^ = rankA'^ 
obtain entries of A-^ by (ITT]) . 

d2(A') = 



/lO -14 
-1 2 

13 -18 

\10 -14 

2. This implies k 



14 
-2 
18 
14 



+ 



+ 



(18) 



1\ 

1 
2 



10 \ 

-1 
13 
10/ 

Ind{A) 



2. We 



10 -14 




10 14 




10 


10 


-1 2 


+ 


13 18 


+ 


10 


10 



2 -2 




2 -1 




18 


13 


-18 18 


+ 


-14 10 


+ 


14 


10 



16 



dii 



3 -14 




3 14 




3 


10 


2 


+ 


4 18 


+ 


3 


10 



and so forth. Continuing in the same way, we get A 



D 



/0.5 0.5 -0.5 0.5 \ 
1.75 2.5 -2.5 1.75 
1.25 1.5 -1.5 1.25 
\0.5 0.5 -0.5 0.5/ 
Now we obtain the Drazin inverse solution x of the system (I18p by the Cramer 
rule (ITB]) . where 



g = aV 



/3 


-4 


4 


3\ 








/10\ 





1 


-1 







2 




-1 


4 


-5 


5 


4 




3 




13 


V3 


-4 


4 






v) 







Thus we have 



Xi 



X2 



^ 8 
1 

Xi = - 



10 


-14 






10 


14 




10 


10 


-1 


2 




f 


13 


18 


+ 


10 


10 


10 


10 






-1 


-2 




-1 


-1 


-1 


-1 


+ 




13 


18 


+ 


10 


10 


10 


10 






2 


-1 




13 


13 


13 


13 


+ 




-18 


13 


+ 


10 


10 


10 


10 






2 


-1 




18 


13 


10 


10 


+ 




-14 


10 


+ 


14 


10 



1 

2' 



1 

2' 
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